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P \ This paper presents a complete analysis of the effects of second order gravi- 

tational perturbations on Cosmic Microwave Background anisotropics, taking 
explicitly into account scalar, vector and tensor modes. We also consider the 
second order perturbations of the metric itself obtaining them, for a universe 
dominated by a collision-less fluid, in the Poisson gauge, by transforming the 
known results in the synchronous gauge. We discuss the resulting second or- 
der anisotropies in the Poisson gauge, and analyse the possible relevance of 
the different terms. We expect that, in the simplest scenarios for structure 
formation, the main effect comes from the gravitational lensing by scalar per- 
q^ , turbations, that is known to give a few percent contribution to the anisotropies 

(— | \ at small angular scales. 
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I. INTRODUCTION 

The increasing number of measurements of Cosmic Microwave Background (CMB) 
anisotropies in the last years and the very ambitious observational programs for the fu- 
ture generation of detectors makes us hope that the angular spectrum of the anisotropies 
will be known with great accuracy within the next decade. This fact has stimulated theoret- 
ical efforts to obtain more precise predictions for the anisotropies produced in the different 
structure formation models, and it is expected that future observations will be very helpful 
in distinguishing among them and in putting constraints on the cosmological parameters. 

Most of these theoretical computations involve numerical or semi-analytic solutions of 
the linearized Boltzmann equation. Non-linear gravitational effects on the anisotropies have 
been computed for some particular processes, such as the gravitational lensing from density 
perturbations ||T]-^J and the Rees-Sciama effect [pHT^] (which is second order in a flat matter 
dominated universe, as the gravitational potential is constant to first order). It has been 
shown that the effect of the gravitational lensing by density perturbations is to smooth the 
so called Doppler or acoustic peaks in the angular spectrum at high £, and it is thus relevant 
for the analysis of the small angle observations [13| . On the other hand, the non-linear Rees- 



Sciama effect is in most cases expected to be much smaller than the first order anisotropies 



14| , |T5H unless early reionization substantially erases the first order anisotropies. 

In a recent paper, Pyne and Carroll |16| have presented a nice framework for a complete 
computation of second and higher order gravitational perturbations of the CMB. Their 
algorithm essentially involves computing the redshift experienced by the photons during 
their travel from the last scattering to the observer in terms of their perturbed geodesies and 
then obtaining the perturbed geodesies up to the required order. The study of second order 
anisotropies is relevant because they can produce a non negligible contribution compared to 
the first order ones, due to the long distances involved in the problem. The reason is that 
several second order terms include integrals of the metric perturbations along the photons 
path, that can enhance small effects as the photons travel from the last scattering surface. 
Moreover, second order effects are also important because they give the primary contribution 
to some statistical measures of the anisotropies that are vanishing for the linear contribution, 
as for example the three-point function of temperature anisotropies [|TT|-|HJ|]. In any case, 
it is important to know the magnitude of the second order effects as they contribute to the 
theoretical error of linear anisotropy calculations. 

In this paper, we apply the formalism proposed by Pyne and Carroll to the computation 
of the full second order anisotropies in the Poisson gauge. We first present a computation of 
the second order anisotropies that generalizes the results of ref . 0] , in that we consider the 
motion of the observer and the emitter, we explicitly include the second order perturbations 
of the metric itself and we take into account scalar, vector and tensor modes. We then 
consider the Poisson gauge, that, in the case of scalar first order perturbations, reduces to 
the longitudinal gauge. We obtain the second order perturbed metric for a dust dominated 
universe in the Poisson gauge explicitly, and then discuss the second order anisotropies for 
this particular case. 

Throughout this paper Greek indices /i, z/, . . . take values from to 3, and Latin ones 
i, j, . . . from 1 to 3. We take, for the metric, signature +2; units are such that c = 1. 



II. TEMPERATURE ANISOTROPIES 

The quantity of interest is the angular variation of the temperature measured by an 
observer. 

We consider a perturbed flat Robertson- Walker space-time and use conformal time 77 
(dr) = dt/a(t), where a(t) is the scale factor of the universe). We can write the line element 

as 

ds 2 a a\ V ) (gVj + gW + Igf) + . . .) d x»dx\ (2.1) 

where the term between brackets is the conformally transformed metric (g^ u ), gffl is the 
background Minkowski metric, and gj£} and gffi are the first and second order perturbations 
respectively. 

Photons travel along null geodesies x fl (X), that we parametrize with A in the conformal 
metric, connecting the observer, at coordinates x@ = (rjo, xq), to the emitting hypersurface, 
that we take at constant conformal time r]£. This hypersurface can be taken as the last 
scattering surface, at redshift zls- At larger redshifts the hydrogen is ionized and Compton 
scattering off electrons (linked to photons by electromagnetic interactions) couples photons 
to baryons. At z^s, hydrogen recombines and photons can start their travel. We assume that 
thermal radiation with temperature Tg(p, d) is emitted by every point with coordinates p; 
in this hypersurface. This temperature depends also on the direction of emission described 
by the vector d, normalized to unity in the background. The different photon paths are 
specified by the direction from which they arrive at O, specified by a vector e normalized to 
unity in the background. This direction fixes the point p and the direction d at emission. 

If the CMB has a black body distribution and the photons suffer a redshift z during 
their travel from the emitter S to the observer O, the emitted frequency uj£ and the ob- 
served one loo are related by uj = uj £ /{\ + z). Since the occupation number per frequency 
mode is conserved, the corresponding photon temperatures are related by T = T £ /(l + z). 
The anisotropies detected by an observer are due to inhomogeneities in the temperature at 
emission and to the different redshift suffered by photons coming from different directions. 
We will compute this quantity up to second order in gravitational perturbations. 

The temperature measured by an observer at O can be written as 

T (x ,e) = ^T £ (p,d), (2.2) 

U) £ 

with co = —g^U^k", where If 1 is the four- velocity of the observer or emitter and k v = dx v jd\ 
is the wave vector of the photon in the conformal metric, tangent to the null geodesic x v [X) 
followed by the photon from the emission to the observation point. In fact, we will propagate 
photons back from the observation point to the emission surface. We thus need to obtain u)£, 
p and d for a given initial set of values Xq, e and 00 q. The resulting quantities are functions 
of the photon path and wave vector, that we expand in series of the metric perturbations 
gj$ and their derivatives as 

a^(A) = x {0 ^(X) + x^(X) + rr (2 ^(A) + . . . , 

P(A) = A; (0) ^(A) + A; (1) ^(A) + A; (2) ^(A) + . . . . (2.3) 



Contrary to the assumptions of ref. 16|| , we are not taking the observer and emitter 



comoving with the total fluid of the universe. In this way we keep track, to the first order, 
of the dipole due to the observer's motion and of the Doppler effect due to the emitter's 
motion, that are otherwise lost. Besides these effects, to the second order, we also take into 
account cross terms involving the velocities and other sources of anisotropy. 
We can expand the four-velocity as 

J7" = - (S$ + v^ + iu< 2 >" + . . .V (2.4) 

This is subject to the normalization condition U^U^ = — 1. 

It is also useful to write the perturbed spatially flat conformal metric as 

Soo = -(l + 2^+^ 2 ) + ...,) , (2.5) 

9oi = 4 1] + U 2) + ■■■, (2-6) 



9* = (l - 20 (1) - (2) ) 5 l3 + x g> + ~ X $ + ..., (2.7) 

whereQ xf — and the functions tp^ r \ z[\ (p^ r \ and x% represent the r-th order pertur- 
bation of the metric. 

The normalization condition for the velocity fixes the time component v^° in terms of 
the lapse perturbation, ip( r >. For the first and second order perturbations we obtain: 

V W° = -V (1) ; (2.8) 

v (2)o = _^(2) + 3 (^(D) 2 + 2z ( t ) v {X)i + vPv™ ■ (2-9) 

In order to obtain the variation in the sky of the observed temperature up to the sec- 
ond order, according to eq. (|2.2|), we need to expand uj and u>£ up to second order in 



gravitational perturbations 

u> = u>V>(l+u,V+u>® + ...), (2.10) 

and also to expand the temperature at emission 

T,(p,d)=Tf (l + r(p,d)). (2.11) 

We will not perform a full expansion of r(p, d), as a calculation of this quantity would be 
beyond the aim of this paper. We will instead assume that it is known for a given model and 
compute the additional effect of gravity along the photons path. We also have to take into 
account that the point p and direction d at emission need to be expanded in the expression 



1 Indices are raised and lowered using 5^ and 6ij, respectively. 



of r(p, d) as p = p^ - 1 + p^ 1 -* + . . ., and d = d^ + d^ + . . .. Performing these expansions in 



eq. (|2.2| ) we obtain [|l| 



T (x ,e) 



(o) 

^O rp(0) 

UJc 



i + (ag> - 4 1} + r 



+ (4 2) - 4 2) + (4 1} ) 2 - 4 1} 4 X) + 4 1} r - 4 1} r 



<9x 



0d* 



(2.12) 



where r and its spatial derivatives have to be evaluated at (p*- - 1 , d^ ^). The first factor gives 
the mean temperature at the observation point Tq = Tg cu /to £ , and the round brackets 
inside the term in square brackets define the first and second order perturbations, ST^ 1 ' and 
8T^ 2 \ that we are looking for. 



To compute them, we will use the same background geodesies as in ref. JL6 

x^ = (A,(A -A)e*), 
jfe(°)'» = (l,- e i ), 

and boundary conditions at the origin 

x^(Xo) = x^(Xo) = 0, 



(2.13) 



(2.14) 



The condition that the wave vector is null fixes the value of k^°(Xo) and k^ 2 '°(Xo). We will 
only need k^°(Xo) explicitly 



k^°(X 



o 



-IpQ Zq 6i 



c(l) , *(!)<*• 



(2-15) 



Using the metric, four-velocity and wave vector expansions we can obtain the quantities 
in the expansion of u 



u® = a~\ 



^) = k^ + ^+vf ) e l + z { pe\ 



u> 



(2) = fc (2)0 + ±^(2) + t z f) e i + L v f) e i 



1 

2' 



I(^(D)2 + \ v <p v w + feCDo^d) _ v w k d)i _ z m k a)i 



- ^zVj - 2<j><-%$ L) e i + xgW^ + 



dk^° 
dX 



AX + (ijf + zfle* + vQe'frM., 



(2.16) 



where AA is the difference in affine parameter between the points where the background and 
first order geodesies intersect the t] = f]s hypersurface, and is given by A A = — x^° + . . .. 
ID that p Wi = x {1)i + x (1)t V and 



It can also be seen 



d^ 



JfcW 

1W\ 



(2.17) 



Finally, we obtain, for the first order temperature anisotropy, 



<.(!) 



W. 



,(!)*« 



,(l)0 .,(!)*< 



,W* 



,(1) 



-0o + -jXo e i e j + v o Ci-ky -vy d-zy d-ipy + r, 



(2.18) 



and, for the second order one, 



ST (2) = ff,(2)0 + 1^(2) + 1 (2) e i + l^) e * _ I ( ^(l) )2 + 1(1)^(1)^ + fc (l)<yi) 

\ Z Z Z Z z 



*V>g>J - 20«,jV + xgVt,™) £ + (4? + 4M* + 



dk^° 



dX 



x 



(1)0 



f + *&V + *&' + r,A { X ™ + xWV) fi + 7 - 



9r 



d«* 



- (kW> + vl 1) e i + z? ) e i + iJ>W-T) e (V 1)0 + v^e* + zfV + ^)\° £ . (2.19) 

The next step is to obtain the null geodesies up to second order; in particular, we need to 
compute the quantities k^°, k^ 1 '^ and x^ 1 ' 1 * to substitute in eqs. ( |2.18| ) and (|2.19| ). This 
problem has been solved for a general perturbed spacetime in ref. |TB| using the geodesic 
expansion introduced by Pyne and Birkinshaw [p0 |. Following their method, we obtain, for 
perturbations around a flat Robertson- Walker background in any gauge, that the first order 
wave vector is given by 



kW{\e) = ^ - 0g } + ixg^ej - 2^ - zf\ % + h(X £ ), 



with 



h(Xs 



x £ 



Xo 



dXA^\ 



where A& = ^ (1) + 0« + zfV - ^xgW, and 



k^(X £ 



>q e + z Xo e j 






with 



AV* 



A f 



dXA {1) '\ 



(2.20) 



(2.21) 



(2.22) 



(2.23) 



For the first order geodesies, we obtain 



x (1)0 (Xe) = (X £ -X o ) 



,(!) AV 



,(*)#« 



y a -<j>o +2 x o iGj 



+ / dX 



-2</> (1) - 4 1} e* + (X £ - X)AW 



(i). 



M l „(i)*i. 



X W* ( A £ ) = ( A , - A ) 2$V + 4T - x^e, 



A f 



tfA 



Xo 



i(i) e < + 2 (i)< _ x im e . + (a £ _ a)^^'*] 

6 



(2.24) 



For the second order, we need only the difference between the wave vector at emission 
and observation 

% ~ ko = TO — Y£ ~ y z £ e i + lj z O e i + ^VO k(D ~ ^TS % 

- (2x^Hf + 2x^ ^y - Z MkV + xWzV'e* + x™ i zf}e?) e + I 2 (X £ ), (2.25) 



with 



h(Xs) = f X£ d\ \\ A W - (4 1} ' - xg^X^ + e^ (1)0 ) 

+ 2fc«°A«' + 2^'A^ + x«°A«" + x^Af] , (2.26) 



where A< 2 ) = ^ (2) + (2) + z\ V - §x£ ; eV 



2 Ar/ 



We can now write the temperature anisotropy in terms of the metric perturbations. 
Replacing eq. ( |2.2(J| ) into ( |2.18| ) we obtain, for the first order, 



8TU = 4 1 ) _ ^g) + v ff% - V W% + r - h(X £ ). (2.27) 

This is a general expression, valid in any gauge, that takes into account scalar, vector and 
tensor perturbations. It includes the effect of intrinsic anisotropies in the last scattering 
surface, dipole due to the observer's motion, Doppler effect from the emitter's motion and 
gravitational redshift of the photons. It is equivalent to the well-known result originally 



obtained by Sachs and Wolfe |21| . It can be seen that the full expression is gauge invariant up 
to a monopole term; the relative contributions from the intrinsic, Doppler and gravitational 
redshift contributions are however gauge dependent. 
Analogously, for the second order, we obtain 

^ = \*f> - \^S + \{^f - \^? - HXe) - v^e^ 

+ (h(X £ ) + 4 1)l e 4 ) (2V#> - $ + lx i o )lJ e l e J - vff% - # - r + vf l e % + h{X e \ 

(i)0j(i)' , (JS)i i JV°J\ f„iM) _ J 1 ) J i ~/\ i .,(!)* f 1 .,^) _o^ 1 ( 1 ).. i J 1 ) J 



+ x^Ag' + (4 ■" + 4 ■' V) {ti - <^ % + rj) £ + v%" {-v% - m 1 * + Xoi,e 



- v§ l ei U$ - (j) { o ] + -Xo^ekej - r - ^ 

+ 4? (-4 1)l + 4? + 2«£V - xTe, - li{Xe)) . (2.28) 

This is also a general expression, that is valid in any gauge and takes into account scalar, 
vector and tensor perturbations. It also includes the effects of the motion of the observer 
and the emitter. In the previous expression we have dropped the terms proportional to v^ % 
as this computation is not aimed at obtaining v l at the emission or observation points, but 
assumes that they are known quantities. 



To proceed further with the computation, we need to know the initial values and the 
evolution of the perturbations. To solve this it is necessary to fix a gauge. There are different 
possibilities: the synchronous gauge (ip^ = z^ r > = 0) turns out to be convenient for many 
calculations and has been widely used for linear anisotropy computations. Another choice 
is the Poisson gauge {zf % = Xij^'^ — 0), recently discussed by Bertschinger [p^] , that in 



the case of scalar perturbations reduces to the longitudinal gauge. The latter gauge, in 

(r) (r) 

which z\ = Xij = 0, has become very popular, because the evolution equations are most 
similar to the Newtonian ones, and thus closest to our classical intuition. All second order 
temperature anisotropy calculations have been performed in this gauge. Since the vector and 
tensor modes are set to zero by hand, the longitudinal gauge should not be used to study 
perturbations beyond the linear regime: this is because in the nonlinear case the scalar, 
vector, and tensor modes are dynamically coupled and vector and tensor modes cannot be 
set to zero arbitrarily. This could be a problem when studying the Rees-Sciama effect that 
explicitly involves non-linearities in the metric perturbations; we will come to this point 
in section IV. We will use the Poisson gauge, which overcomes the above limitation of the 
longitudinal gauge, while keeping all its advantages in terms of physical interpretation of 
the results. 



III. SECOND ORDER PERTURBATIONS IN GENERAL RELATIVITY 

We consider the gravitational instability of irrotational collision-less matter in a flat 
Robert son- Walker background up to second order. Different approaches to this problem 
have been proposed. The first solution of the second order relativistic equations has been 
obtained, in the synchronous gauge, in a pioneering work by Tomita p3fl . Matarrese, Pan- 



tano and Saez p4|,[25| obtained the leading order terms of the expansion, using a different 



method, based on the so-called fluid-flow approach. Salopek, Stewart and Croudace |26 
used a gradient expansion technique to obtain second order metric perturbations; an intrin- 
sic limitation of their method is, however, that non-local terms, such as the non-linear tensor 
modes, are lost. Russ et al. |27j recently rederived the metric perturbations to second order 
in the synchronous gauge, using a tetrad formalism. We are interested here in obtaining the 
solution in the Poisson gauge. Instead of perturbing the Einstein equations in this gauge 
and then solving them, we will transform the solution known in the synchronous gauge to 
the Poisson one, using the second order gauge transformation recently developed in ref. EH 
(for more details see ref. p5fl). 

Up to this point we have been completely general in the inclusion of scalar, vector and 
tensor modes. Now, in order to give a more quantitative insight on the relevance of the 
different contributions, we will make some restrictions. First, we will neglect vector modes 
at the linear order. The fact that they have decreasing amplitude and that they are not 
generated in inflationary theories, makes us expect that their role will not be relevant, at 
least if inflation was the mechanism for primordial fluctuation generation. We will however 
keep track of the second order vector modes generated by the coupling with scalar modes. 
Second, we will neglect the effect of linear tensor modes as sources for second order metric 
perturbations: because of the graviton free-streaming inside the horizon this is a very rea- 
sonable approximation. We will then only consider the second order perturbations generated 
by linear scalar perturbations. These are expected to give the dominant contribution, as 



vector and tensor modes have decaying amplitude. The solutions given in refs. f23,27fl, and 



hence those obtained in this section apply to this restricted case. For the computation of 
the CMB anisotropics we will however keep the contribution of the linear tensor modes in 
ST^ and ST^ everywhere, except as source for gj?J. The contribution of tensor modes can 
in fact be comparable to the scalar contribution to ST^ at large scales in many inflationary 
models [fllj . 

The first order solution to the perturbed Einstein equations in the synchronous gauge is 
given by (see, e.g., ref. 



4" = 4" = o, 
4"^4vv 

7 2 ( 



,(D _ ^ / 1, rr2\ , T(i) 



X K & = -\ ^-i««VV +xi w , (3-1) 



T (!) W 



where (p = <p(x) is the initial peculiar gravitational potential. Xij ls t ne tensor (transverse 
and traceless) contribution that can be written as 

X J (1) (x,77) = ^|ci 3 kexp(zk-x) X «(k,r / ) e -.(k), (3.2) 

where e?-(k) is the polarization tensor, with a ranging over the polarization components 
+, x, and xP (k, if) is the amplitude. Its time evolution during the matter dominated era 
can be represented as 

£\K v )^A(k)a a (k)(^^y (3.3) 

where a CT (k) is a zero mean random variable with auto-correlation function (a (r (k)a -'(k / )) = 
(27r) 3 fc~ 3 5 3 (k + k')5 a a>- The spectrum of the gravitational wave background depends on the 
processes by which it was generated, and for example in most inflationary models, A(k) is 
nearly scale invariant and proportional to the Hubble constant during inflation. 
The second order perturbations are given by [p3|, p7| , p^ 



4 2 > = 4« = o, 



l2> ~ h. (-f^-"' + (vV)2 ) + t/ (^ + r v ' v 



's 

) _ 

3 

5 o / „ , . „ u „ . 4 



(2) _ V fm„A„ 10,„ V72,„ , >im2,^2jc 19 



X% = 7^ (/ 9 W,jy - 12<^;V> + 4(W)^« - y ^VM. 



+ ^ 2 f -6<P,i<Pj ~ bpH>,n + 2<P' k( P,k5ij + -<pV 2 <p5 i:J ) + tt^, (3.4) 

where the traceless and transverse contribution rcsij satisfies the inhomogeneous wave equa- 
tion 

An + -An - V 2 vr^. = -|^V 2 ^-, (3.5) 

9 



with 

Sij = V 2 * <% + *o« + 2 (^,y VV - <P<ik<Pj) , (3-6) 

where 

V 2 *o-^((VV) 2 -W ifc ). (3.7) 

This equation can be solved using the Green method; we obtain for irsij that 

ri 4 4 

TTSr,- = —Sij + -77 2 Ty(x) + fly (x, 7/), (3.8) 

where V 2 Tjj = SV,- and the remaining piece ffy(x, 77), accounting for a term that is constant 
in time and another one that oscillates with decreasing amplitude, can be written as 

~ / \ 1 f jSi /-i \40„ ,, s /l cos( kr]) sin( kr])\ 

-«(*, ») = p^i / A«p(4 ■ x)^(k) (3 + ^ - -jj^J , (3-9) 

with -S'y(k) = /<i 3 xexp(— ik • x)S'ij(x). 

The gauge transformation is determined to each order by a four-vector £( r )^ that we split 
as £( r >° = c*M and £( r >* = ^W + rfW*, with d^ = 0. In ref. |2B| the vectors f ^ and 
£(2)/^ describing the gauge transformation from the synchronous to the Poisson gauge, have 
been explicitly obtained in terms of the synchronous metric perturbations g^ and g$L- 
Using the metric perturbations in the synchronous gauge presented above, we can write the 
first order gauge transformation as 



(i) V 
3 



a" ' = -ip, 



P (1) = £<p, (3.10) 



and dP~> 1 = 0, in the absence of vector modes in the initial conditions. For the second order, 
we obtain 

a (2) = -^V 3 ^o + V (y^ 2 + 46 



P w = V 4 ( ^'V,< ~ ^0) + -f [^ + 66o ) , (3.11) 



where V 2 6 = ^0 — \f ,l f % an d 



4 _ 4 



V 2 df = rf [--^-VV + ~^,ij ~ ~*o,;J • (3.12) 

We can now compute the metric perturbations in the Poisson gauge using the transfor- 
mation rules of ref. |28 . For the first order, we obtain 

4 ] = 4 1} = v, 

X% = xZ 0) . (3-13) 

10 



These equations show the well-known result for scalar perturbations in the longitudinal 
gauge and the gauge invariance for tensor modes at the linear level. For the second order, 
we obtain 



# 


= V 2 QvV,« - ^*o) + fv 2 + 126 


# 


= V 2 (~p'V,i - ^o) + 3V 2 - 86 , 


v 2 4 2) * 


= -^(^VV-^VjH-2^)> 


(2) 

Xpij 


= ftij- 



(3.14) 

Note that the resulting expressions for ipp and cj)p can be recovered, except for the sub- 
leading time-independent terms, by taking the weak-field limit of Einstein's theory (e.g. 
ref. [p2| |) and then expanding in powers of the perturbation amplitude; this is basically the 



method employed in previous second order computations of the Rees-Sciama effect. Also 
interesting is the way in which the tensor modes appear in this gauge: the transformation 
from the synchronous to the Poisson gauge has in fact dropped the Newtonian and post- 
Newtonian contributions, whose physical interpretation in terms of gravitational waves is 
highly non-trivial (see the discussion in ref. |£J); what remains is a wave-like piece plus a 
constant term which has no effects on ST^ 2 \ 



IV. POISSON GAUGE ANISOTROPIES 

Let us start by discussing the first order anisotropies that are described by eq. ( [2.271) . 
The first term ip £ represents the contribution from the gravitational redshift of the photons 
due to the difference in gravitational potential between the emission and observation points. 
iPq only contributes to the monopole and can be neglected. The term v e, is the dipole 
due to the motion of the observer. The term v £ et accounts for the Doppler effect due to the 
velocity of the photon-baryon fluid at recombination and contributes to the acoustic peaks. 
The term r describes the intrinsic anisotropies in the photon temperature and is highly model 
dependent. For example, for adiabatic perturbations, in which all the components (baryons, 
photons, dark matter) have a constant number density ratio, the photon energy density, and 
thus the temperature, varies proportionally to the potential fluctuations (at scales larger 
than the Jeans length). It can be seen that in this case r = \ — \s = \- £l - = — §<%■ • It is 
the combination of this term and the first one that gives the standard result for adiabatic 
perturbations at large angular scales, 5T = |y?. At small scales, r gives the main contribution 

to the acoustic peaks. We have not intended here to include a computation of r and v £ , 
that would involve solving the linearized transport equation for the photons, that is coupled 
to the fluid evolution equations for the cold dark matter component and the baryons, the 
Boltzmann equation for the neutrino distribution and the Einstein equations for the metric 
perturbations. This problem has been treated and solved numerically by several authors 
(see, e.g., P^j- fJOfl ). We assume that r and v £ are known for a given model, and compute 



the additional anisotropy generated by the metric perturbations along the photon path up 
to second order. 

11 



Finally, the contribution to 5T^ from the last term, Ii(\g) (called integrated Sachs- 



Wolfe effect and given by eq. (|2.21|) ), represents the additional gravitational redshift due 
to the time variation of the metric during the photon travel. As in the linear regime the 
scalar potentials and if) are constant in time for a flat matter dominated universe, their 
contribution vanishes (they will however give a non- vanishing contribution in the non-linear 
regime). The contribution of tensor perturbations to the temperature anisotropies arises 
exclusively from Ji(Ag) at the linear order. In many inflationary models, in which besides 
the usual scalar perturbations also a background of gravitational waves is produced, their 
contribution to the CMB anisotropies can be comparable to that of scalar perturbations at 
large scales. The contribution from the observer's motion term is of order 10 -3 , while the 
remaining part contributes for an order 10 -5 . 

We can now discuss the second order anisotropies that are given by eq. ( |2.28 ). The first 



Then, there is the term (x £ + x £ e 3 ) (ipj — v\Je l + t\j J , that can be split into a 



(2) 

term, given by ip £ , represents the gravitational redshift of the photons due to the second 
order metric perturbations, and is much smaller than its first order equivalent. Then, there 
are several terms involving products of two of the terms contributing to ST^ ; these are all 
very small compared to ST^ l > (at least three orders of magnitude smaller) and can safely 
be neglected. Also the term (dr/dd z )d^ z is the product of two small quantities and can be 
neglected. 

,{i)j , „.(i)o„^ fj,W „,(i)„i i J 1 ) 

piece proportional to 

= (A, - Xo) (-xT 3k e k + xT ik e k e^) + f* dX ( X ^ k e k - X T{1)lk e k e^) 

- £' dX(Xs - A) (V " 2^,*eV - ±x£ (1)J eV + ^xgJWeV) , (4.1) 

and another piece proportional to 

(xj 1)j + 4V) £ = e j (xi 1] e l + x^^e 

= -e> £' d\ (ty - IxTe'e 1 ^ . (4.2) 

The first piece describes the effect of the gravitational lensing on the photons as they travel 
from the last scattering surface to the observer. The transverse displacement x]_ includes 
the usual contribution from scalar perturbations ((f), that has been considered in some 
previous studies fT3|,|6| and has an observable effect on small angular scales, and a new 
contribution due to the gravitational wave background (xij ) acting as a source of lensing. 
The second piece is due to the time delay effect of the lenses that changes the spacelike 
distance to the intersection of the photon path with the last scattering surface. The scalar 
part of this term is expected to be suppressed with respect to the gravitational lensing term 
due to the spatial derivative of ip that appears in eq. (f4.1|). The gravitational waves part is 
probably of the same order of magnitude as its gravitational lensing counterpart. The term 
x £ A £ is similar in form to the gravitational lensing and time delay terms: it arises due 
to the difference in affme parameter along the perturbed and background geodesies. As it 

12 



involves time derivatives of the metric perturbations, the contribution from scalar terms will 
be small, but the tensor contribution is expected to be larger. 

The next term to consider couples the velocity of the photon-baryon fluid with the per- 
turbation to the photon wave vector at emission and is given by v s ^I\(\s). Comparing it 
with the gravitational lensing contribution, we expect some reduction because the Doppler 
contribution to the first order anisotropies is smaller than the other first order terms (al- 
though of the same order of magnitude) and some enlargement because of the factor (Ag — A) 
of difference between the I\{\s) and x ± expressions. A more careful quantitative estimate, 
which would require a choice of the particular structure formation model of interest, is 
beyond the aim of this paper. 

Finally, we have the term ^(Ag), that is given by eq. (|2.26|) and is an integral of several 
terms. The first one accounts for the Rees-Sciama effect, given by 

STns = \ £ dX (>' + ^' + zfj - i X Sf eV) . (4.3) 

We can use the second order perturbations of the metric obtained in section III to compute 
it. The contribution from the scalar perturbations, ip^ and (f)^ 2 \ is given by 

rvs /l • 20 \ 

6T RS = J d m {-<p>Yi ~ ^ *oJ , (4.4) 

where the terms inside the brackets have to be evaluated along the background geodesic 
parametrized by A = t]. This piece coincides with that considered in some previous studies 



of the Rees-Sciama effect [H- 10]. The resulting anisotropies turn out to be between one 



and two orders of magnitude smaller that the first order ones. The contribution from the 
vector and tensor modes can be obtained by substituting z\ and xij from eq. ( J3.14 ) into 



eq. (fPf ). Let us estimate their magnitudes compared to that of the scalar piece. The 
integrand for the vector piece is z\ e % ~ kip 2 , while the scalar one is ip™ ~ r]k 2 Lp 2 . Thus, 
z\ e l ~ ?/> /{f]k) ~ ?p( 2 > (aH/k) and the vector contribution is suppressed with respect 
to the scalar one, as the wavelengths of interest are smaller than the Hubble radius. This 
estimate is similar to the one obtained in ref. |14|]. In the same way, the integrand for the 
tensor piece is xij e%e ? ~ k(p 2 /(kr]) 2 ~ ip^ /(kr]) 3 ~ ip^ (aH/k) 3 . Hence, also the tensor 
contribution is much suppressed with respect to the scalar one. 

The integrand of the second term contributing to -^(Ag) is Xij e^(k^ 1 + e l k^°), it 
represents a correction to the anisotropies generated by linear gravitational waves, due to 
the perturbation of the photon wave vector. The piece containing k^- 1 ' is expected to be 
smaller than the other one; the largest contribution can arise from the term 

X£ rfAxJ^V I" d\'AU>\ 

■o JXo 

The last four terms in ^(Ag) will have a small contribution coming from scalar per- 
turbations as they involve time derivatives of the potentials ip and <fi that are constant at 
linear order. The contribution coming from gravitational waves is expected to be larger, in 
particular the last two ones 

J x yx( X U%W'eV+ X M%M'e^). 

These can be interpreted as the gravitational lensing and time delay effects acting on the 
anisotropies generated by the linear gravitational wave background. 
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V. CONCLUSIONS 

We have computed the anisotropies in the CMB radiation up to second order perturba- 
tions in the metric around a flat Robertson- Walker spacetime. This calculation generalizes 



the results of ref. [16] in that we have taken into account the velocity of the emitter and the 
observer, we have considered scalar, vector and tensor perturbations and we have explicitly 
included the second order perturbations of the metric. We have obtained these second order 
metric perturbations for a universe filled with a collision-less fluid in the Poisson gauge, by 
performing a second order gauge transformation of the synchronous gauge solutions, that 
have already been studied in some detail in the literature. 

Using these results, we have discussed the relevance of the second order contributions 
to the anisotropies in the Poisson gauge. The most relevant expected contribution is due 
to the gravitational lensing of photons due to density perturbations, that has already been 
the subject of several studies. We have shown that also a gravitational wave background 
acts as a source of lensing for the CMB photons. This effect is much smaller that the scalar 
one for a gravitational wave background with spectral index ny = as generated during 
an inflationary period. Other sources of gravitational waves with more power than the 
inflationary ones at small scales may give a larger contribution through this effect. Other 
contributions include the time delay effect of scalar and tensor lensing, a coupling of the 
velocity at emission with the perturbed photon wave vector and a second order perturbation 
to the integrated Sachs- Wolfe piece. This term includes the well-known Rees-Sciama effect, 
that has been widely studied for the time variation of the scalar gravitational potential. 
Using the second order perturbed metric in the Poisson gauge obtained in section III, we 
have shown that the additional contributions to the anisotropies arising from the vector 
and tensor modes induced by linear scalar perturbations are expected to be suppressed with 
respect to the scalar one. We have also pointed out the existence of two more terms that 
are corrections to the anisotropies generated by the linear gravitational wave background, 
due to the perturbation of the photon wave vector and to the lensing and time delay effects 
on gravitational wave anisotropies; these can give a relevant contribution to the integrated 
term. These contributions deserve a more detailed quantitative analysis. 

Although in the light of the present analysis we do not expect that the second order 
gravitational effects will give a major contribution to the anisotropies at any scale, it is 
interesting to know if they could be detected by the planned high accuracy satellite obser- 
vations. The gravitational lensing by scalar perturbations is known to give a few percent 
effect in some structure formation models and thus will be relevant if the expected 1% sen- 
sitivity is achieved. The amplitude of the second order terms is also important because 
they contribute to the theoretical error of the anisotropy computations that will be used to 
determine the cosmological parameters from the measured multipoles. 
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